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Об одном подходе к тестированию новьїтх методов решения 
нестационарной задачи теплопроводности 


У роботі запропонований метод побудови точних розв'язків нестаціонарної задачі теплопровідності в 
областях, складених з прямокутників. Пропонується аналіз результатів обчислювального експерименту, 
проведеного авторами за допомогою створеної ними програми при розв'язанні прикладу в областях, 
складених з прямокутників. Розглянуто приклад області, що має форму швелера. 

Ключові слова: нестаціонарна задача теплопровідності, інтерлінація, 

метод скінченних елементів, початково-крайова задача. 


Іо Фіз рарег, Фе пефоа Їог сопзігисйпе ехасі 50/ийоп5 ої ипзісаду Пеає сопаисйоп ргобіет Їог агеа5, 
уурісі аге пзаде ої гесіап5іез, 15 ргорозей. ТПе апаїузіз ої Ше гезиіз ої сотриїанопа! ехрегітепі 15 оЙегед. 
Ії мав сопдйсіва бу Ше ашіогя мій Бе Пре!р ої бо ргоєбгат5 Їог 50Іміпе Ше ехатріе їп агеа8, мрісі аге 
таде ої гесіап5іез. ТРре ехатріе ої агеа, урісП баз а 5раре ої сраппе! Баг5, 15 соп5ідегей. 

Кеу ууогд8: ипзсеаду Беаї сопдисйоп ргобіета, іпіегійпайоп, Ппіїе еЇетепія пешод, 

Фе іпічаї-боппдату уаше ргобіет. 


В работе предложен метод построения точньїжх решений нестационарной задачи теплопроводности в 
областях, составленньх из прямоугольников. Предлагаєтся анализ результатов вьтчислительного 
зксперимента, проведенного авторами с помощью созданной ими программьт при решений примера в 
областях, составленньх из прямоугольников. Рассмотрен пример области, имеющей форму швелера. 
Ключевьге слова: нестационарная задача теплопроводности, интерлинация, 

метод конечньх злементов, начально-краєвая задача. 


Вступ 


У праці (1 | запропонований метод ЛІДР розв'язання диференціальних рівнянь з 
частинними похідними (метод зведення до системи лінійних інтегро-диференціальних 
рівнянь). У статтях |2|, |3| досліджувалися деякі аспекти чисельної реалізації методу, 
запропонованого в статті |1). У даній роботі пропонується аналіз результатів обчис- 
лювального експерименту, проведеного авторами за допомогою створеної ними про- 
грами при розв'язанні прикладу в областях, складених з прямокутників. У цьому 
прикладі розв'язувалась тестова крайова задача, в якій був відомий точний розв'язок, 
побудований з використанням методики, викладеної в роботі |4|, відповідним чином 
перенесеної на випадок нестаціонарної задачі теплопровідності. 
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У статті |2| розглянуто більш детально структура наближеного розв'язку, коли 
областю є прямокутник. У даній роботі досліджуються його можливості для випадку 
областей, складених з прямокутників. Зокрема, ідея цього методу використовується 
для побудови точного розв'язку нестаціонарної задачі теплопровідності, що належить 
класу СР"! ( р х (0, чо )) . Розглянуто приклад для швелера. 


Основні твердження методу 
Наближений розв'язок задачі: 


/ 


Ци(х, у пу м - (ад (х, уди!), (а, (х, ушу), - ГО, у, І), 


(хууєОск,, 150, (1) 
и(х, у,1)-0, (х,у)є дО, (2) 
и(х, у, 0) -Ф.(х, у, (с у)є0, (3) 
Фо. У 70, (4) 


де О - область, складена з прямокутників - шукається у вигляді: 


шобубі У У Сі С) Апх-Ю пу - у) 


ке! узі 


п-1 по-1 п-і п-і 


"У У С.) Кп'х - і) Нпу- 0 - У У С, нах -Ю Н(пу- 0), 


Феї ізі Кар о Фа 
вякій С б (1) - невідомі функції, які задовольняють початкові умови 


Ко)! й Й РОДУ 
Си, (0)-ф б 2 , Су. (0) - ра 2) є з Мп п (0)-Ф р 2) ; 
поп п' п п 


Для їх знаходження розв'язується задача Коші для системи звичайних дифе- 
ренціальних рівнянь відносно С - (1): 


0 
ре РА ОУН ОО (7 Фа. 0) 


(арт аїр У вбкщ к.е-п-1), 


а жо о Нилюсїоз (7 Фа) 


(й -лл? жо к,Ф-1п-1), 


0 о 
р» РР ПУ (говні (К, Р -1,п- 1), 


де використовуються позначення: 


ди, бу 
-,|Би(х, у,1)у(, у,2)4хау, (и, у1( --7--|4хау. 
Цехи ь змо і -Щаб сен ча и й 
Метод оснований на заміні формул інтерлінації |5| за змінними х, у відповід- 


ними формулами інтерполяції з точністю, яку має метод ЛІДР. 
Згідно із загальною теорією методу скінченних елементів (б, с. 290|, якщо по- 


п : й 
значимо через 5" простір МСЕ степеня (к - 1) за кожною змінною, то похибка апро- 


ксимації методом Гальоркіна задовольняє нерівностям: 
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(0) «(о Ра че 


(52 
о, -ежо9д и (є ЇЇ «| 
0 


де п - максимальний крок и С - деяка стала, 


зоні но рено (2) (2) парк 


Для цієї формули справедлива така оцінка похибки в енергетичній нормі: 
і і 2 р 2(к-і 2 
«(и-8'и-8")«С р Чаї, 


де К - степінь скінченних елементів, 


уч ди ду посух бе: 
Ох Ох 


| ( (1) се з 
о 


Якщо ж використовуємо метод скінченних елементів, побудований на основі 
методу інтерлінації (ІМСЕ), то формула (5) для похибки зміниться і буде мати вигляд: 


ї йа| , (6) 


«сі б не 


(и (г) сей (2 «Си Г; (2 фей и | «Деко и 
0 


Постановка задачі 
У даній роботі область )) будемо вважати розміщеною в квадраті Д - |0,1|х|0, 1|, 


складеною з об'єднання прямокутників Й, - Гах, Їх су» А |; К -1,М., зі сторонами, 


паралельними осям координат. Проведемо розбиття на прямокутні елементи за допомо- 


гою прямих хаех,, Кк -0,т, уз у, 2 -0,п так, щоб Уа,, Р,, с,, д, серед точок роз- 
биття існували такі хи» Ук» Чо Ук ЩО Ру» В Ху» С З Уціу» ЧУ 


Треба знайти наближений розв'язок и (98 у,г) початково-крайової задачі 


один (зм); (зу)єр,і»о, (1) 
и(х,у,0)29,(х, У), (8) 
«(хг -0 (9) 


методом скінченних елементів, викладеним в статті |1|, і порівняти його з наближеним 
розв'язком, знайденим МСЕ, в класичній формі. У роботі |3| досліджувався метод, 
викладений в |1|, для прямокутної області, в даній роботі наведемо результати для 
області, що має форму швелера. 


Побудова точних розв'язків тестових задач 
в областях, складених з прямокутників 


Приклад. Розглянемо задачу нестаціонарної теплопровідності в області Р, яка 
має форму швелера (рис. 1). р П АП,,де П, «|-а,а|х(0,Ь|, ЛП, «|-с, се |4,Ь1. 
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Рисунок 1 - Область Р у формі швелера 
Розіб'ємо область |) на прямокутні елементи прямими хех,, К- 0, т; 
узу, Ге0,п; 


зазмєкєесх, т СоХ 


зі чаю, ЗО но 


ть чі 
О-укуко3я Уа -ах« Уа зно Уа, еф, 


наржежения розв'язок поставленої задачі шукаємо у ВиРуда 


т-1 


ім у,1) мУХо 0 (зн У Ус Он роз Ус,(г )Н, (У), 10) 


ат (епі 
пай но 120, 
о. од 
ХХ 
зон хх охох,, 
Хх.-Х 
х кл 
й, оре 
ХХ 
точне, хАхАХх 
Ху 
0 а 


Аналогічно визначаються функції Н, (у). 
- КК єкт ІсИс по ІСЕЄ т п рєИХл-Ї тк тої пс а- 1) 
Позначимо: 52 5,0 3,9 3,, де 3,2 (6 )рскст-Г 1хЙхл,-ії), 
«ПРА ех пу пєбхл-1), З, є Б )|в, іх сто п хЙКа-Щ). 


Функції С,, (1) (о, / с 5) знаходимо, мінімізуючи вираз (п) з використанням 


методу локального потенціалу: 


«(а умах 


Це приводить до системи лінійних диференціальних рівнянь відносно С,, (1) 


Запишемо (1) у вигляді: 


и(х,у,г)з у3 Си (г) в, (х)Н, (у). (11) 


(к.е З 
Тоді для визначення чі )отримаємо че 


доти си (0 Р ии Є й г) й їі 1» 0, (р,д)є б (12) 
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де 


ас З (28, Оп (Н, дахау, 
имя дн Он уні б Ок ан аку 


Й (1) ьо Н, (У) 7( о у,г)дхау. 


Проведемо сумування по кожній множині 45, 25,, 5, окремо. Отримаємо замість 
(12) систему о У зав 


25 ана Сн 2 ан 2, «о е 


(о Гуз, (2,0)еб, (к.е 


(13) 
з 2. Брак Си (2)- ), Баки Си (2)- 53 Бані Си (уч бі (2), 20 
(к.0єЗ, (к.0ез, (к ОЗ 
Си (0) -Ф(х, у,) : 
Систему (13) запишемо у матричному вигляді: 
певно » 0, (14) 


р 


вільних членів /,, (г) на вектор С (г), матрицю С,, (г) на вектор с(!) згідно з опи- 


на матриці А В, та матриці 


шляхом заміни ЧОТИривимірних масивів а рут цу 


раки? "ракі 


саним нижче алгоритмом. 
Тут при формуванні результуючої матриці А, що відповідає написаним вище 
трьом матрицям, поводимося наступним чином: 


А З А» Іхш «(т -1)(п, -1), Ісу «(т -1)(п, -1), 
-А2,у» 1Хр.Х (т, - 1)(п-1- п), су (пт, -1(п-1-п,), 
з АЗ 


(т-1Уь 1 (топ х ку 
У НУ 1)( пу -1)-К(ть -1)(п-1-пу)чи, (т-1)( пу -1) (т; -1)(п-1-пу )ям 
Іси б Ісу «(т, -1-т,(п-1-п,) 
Аналогічно формуємо матрицю В. 
Зауваження. Якщо а -р, то в результаті отримаємо систему Ас-- Вс - Ї, яка 
досліджувалася в |3| для прямокутної області. 


ум? 


Побудова точного розв'язку тестового прикладу 
з використанням сплайн-інтерлінації функцій 
Побудуємо точний розв'язок задачі для тестового прикладу. Необхідно, щоб 


точний розв'язок належав до класу функцій, які мають неперервну похідну по / і не- 
и" 2, 2 
перервні другі похідні ш/,; и/, - Для цього будемо шукати точний розв'язок тестового 


прикладу у вигляді 5 різних формул у 5 різних підобластях Д,, Д),, В, Ду, ДО, 
розбиття (рис. 2). 
Будемо вважати, що на лініях у-0, у-а, у- Р точний розв'язок визначається 
наступними формулами: 
д'ш 
ду" 


5 


ди д'и 
- Фо: (51), 5. Фі), аг фо), уза) 2. 


ує0 узі уєб 
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цецао (1) 
чаша») цу З Ча1(»0) пацо(»,) 
чаші) Фах) 


шу З а (1) 


Рисунок 2 - Розбиття області Р) на підобласті та граничні умови 
на частинах границі області 0, де 0-0, ор, 2РруоДрДус-б, 


Точний розв'язок налініях х--а, хое-с, хаос, х- а визначається формулами: 


Ой! Ой! Ой! 
| Зо (Ум), | зм (ум); | зу» (Ум); 
д'ш 
1 М мг, (у), 5 20,2 
При цьому 
фо(зі) зтачхазе фи(ьг), за«хач-є 
Фіо(г) з 0,-сах«с й Фаїла)в 0,-сахас , 
ФС), сах«а Фи(ог), сахх«а 
фо(5ї), -аєхач-є 
фо дуве 0, -сах«с й 
Фоїні,сєхей 
У о(71) 0«ух«д Уо()1) 0«ух«а 
,1)с і ; ,1)- і з 
нм) і Ум) 0а4«у«Б 
У (1), суха Ум), «ух 
у - 5 7 ,Ї -- і й 
"0 на "0 о за 
у КЕ о«ух«а ТД 9 Р Ор 
:(Ї)з Е зПДузії- , 
шо я на Ом) 0, д«ухь 


де Ф,, (х,1), ко 0,2; М (»,1), «0,3; 045,р42 - довільні функції, які задо- 
вольняють умови С.М. Нікольського (5| у кутових точках. 
Зокрема, в точці (-а,0)мають вигляд: 


ФС аг) 2(0,), Ох ура, 
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що є наслідком незалежності мішаних похідних від порядку диференціювання: 


д""ш д""ш 
Рарерчирі НН ан Х0ювред. 
дх" ду ба ду'дх сао) 
Згідно з твердженням праці |7), у кутових точках області ) точний розв'язок 
зеоьій б | -. ЛФ 
крайової задачі для диференціального рівняння Лапласа має особливість г? "5іп --- 
а 


бу г(х, у) - відстань від точки (х, у)є 0 до кутової точки, ( - вхідний кут, ф 
полярний кут). 

Враховуючи однорідність граничних умов, в загальному випадку частинні 
похідні першого порядку в кутових точках (-а,0), (-а,Б), (-с, Б), (с,Б), (а, 5), 
(а, 0) повинні дорівнювати нулю. Мішані похідні в цих точках не дорівнюють нулю. 


Частинні похідні другого порядку у цих точках теж дорівнюють нулю. Мішані 
похідні другого порядку у цих точках не дорівнюють нулю. 


У точках зі вхідними кутами, тобто в точках (-с,Ф)і (с,4), розв'язок и (х, у,г) 
дорівнює нулю, в загальному випадку частинні похідні першого та другого порядку і 
мішані похідні (з, у) (0х ра 42) не дорівнюють нулю. Але у частинному 
випадку вони (або деякі похідні) можуть дорівнювати нулю. 

Точний розв'язок и (х, у, г) будемо шукати у вигляді: 

ХУ), 12, у)є ), 

ін(х У»), (х,у)є О,, 

(ху) з 3йу(х У, 0), (х,у)Є Ду » 

У ЦЕ Б А М У о 

ТР ОО РО РИ Я й 

де функції и, ц-1,П,ПІ,Г/,/ задовольняють умови на границях між сусід- 


німи областями ),, що забезпечують виконання властивості и є Зо (в). 


Зокрема, 


К 
2 2 1 2 
Во ООН 


кеї ред «0 де0 

При побудові тестового прикладу в даній роботі будемо використовувати лише 
два підходи (для побудови допоміжних функцій). 

Підхід 1. 51, (У) будуємо у вигляді поліномів 5-го степеня, якщо вимагаємо, 
щоб у формулі інтерлінації використовувались і сліди функцій, і сліди перших та дру- 
гих похідних (в цьому випадку по змінній у маємо всього 3 вузлових точки у, 20, 
у 24, у, « Б). Якщо у формулі Ерміта використовуються сліди функцій і перших 
похідних по у, то кожна така допоміжна функція буде поліномом 5-го степеня (6 ко- 
ефіцієнтів можуть бути однозначно знайдені з 6 інтерполяційних умов). Якщо ж ви- 
користовуємо формули, які при уз у,, у з Уу,, у з у, мають задані сліди та задані 
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сліди першої і другої похідних, то потрібно використовувати для побудови допоміжних 
функцій 51, ( у) поліноми степеня 8, які мають 9 коефіцієнтів. 


Аналогічно по змінній х. Якщо в кожній точці інтерполюємо лише функцію і 
похідну, то потрібно використовувати поліноми 7-го степеня, оскільки в 4 точках буде 
задано 8 умов. Якщо ж у кожній точці інтерполяції по х задана функція та перша 1 
друга похідна, то на 4 точках по х нам потрібно задовольнити 12 умов. Таку кількість 
умов можна задовольнити поліномом 1 1-го степеня. 

Підхід 2. Якщо використовуємо як допоміжні функції сплайни, то при умові, 


В МР 7 У ДР Че 
що нам необхідно, щоб ТОЧНИЙ розв'язок належав класу С (0) , Достатньо взяти сплайни 


3-го степеня, оскільки в цьому випадку на кожному інтервалі (х,,х,.) в точці х, та 


в точці х,,, можемо відновити значення функції та її перших похідних (4 умови). 
. УЛ Р . 
Якщо хочемо, щоб тестова функція належала класу С (0) , потрібно викори- 


стовувати сплайни 5-го степеня, оскільки в точці х, та в точці х, , будуть задані зна- 


чення слідів функції, першої та другої похідних (6 умов). Тобто на цьому інтервалі 
можемо скористатися поліномом 5-го степеня, який має 6 коефіцієнтів. 


Функції 21 РИ (3), 52 пр (є), Ко 0,2, і- 0,3, ре 0,2 є відповідними базисними 


поліномами Ерміта, їх похідні відповідного степеня повинні задовольняти визначені 
умови у вузлових точках. При визначенні цих функцій використовується позначення: 


бе 5 бут. 


Зокрема, 510 ( у) - поліном 3-го степеня з властивостями: 


451 6 (У "51 5 (У 
СІРА БУ Б 5 шо (У), 70; і й ен -0. 
у-0 ує0 
ор 1 у у? 
51 не -4 учнанньвея ши - м ньсеАни 
обучочеуутеі тона 


Припустимо, що точний розв'язок задачі має вигляд: 
У 
и(х,у зе" й(х, у). 
Права частина точного розв'язку задачі (7) - (9): 


Ї (- у.) зей" (2 й(х, у) з ді (х, у) 


Графік функції / фо у,г) при 1-0, 2. - 2, яка визначає щільність розподілу точ- 


кових джерел тепла у пластині, має вигляд (рис. 3): 
". 20 


Рисунок 3 - Графік функції / (, у,г) при 1-0, 2-2 
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Про один підхід до тестування нових методів... 


Проведений обчислювальний експеримент показав, що відхилення набли- 
женого розв'язку початково-крайової задачі з вказаною правою частиною 


(є у, г) ай (-х й, у) -Ай(х, у) має в рівномірній нормі похибку є - о(п?) , де 


п - крок розбиття, згідно з теоретичними твердженнями, що задовольняє спів- 
відношенню (6). 


Висновки 


Таким чином, запропонований у роботі метод побудови точних розв'язків неста- 
ціонарної задачі теплопровідності в областях, складених з прямокутників, дозволяє 


використовувати ці розв'язки для випадків, коли и (х, у, г) ебо- (р х (0, ко). 
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О.М. Іупууп, І..9. Т.обапота, С.У. Хаїугппа 
Оп те Арргоасп іо Тезппо ої) Меу Метоа5 /ог 5о0Їуїпо 
те Цпзтеаду Неаї Сопаиспоп Ргофбіет 


Трі8 рарег рге5епіз ап апаЇузі8 ої Ше гебиіїз ої питегіса! ехрегітепіз сопдйсіва 
Бу Фе апШог8 ипзіпя, Ше ргоєгат, урісі ШФеу сгеаїе іо 50Їме ргобіет5 їп агеаз5 соптпро5ей 
ої гескапеїе5. Па (Ві5 ехатріе, Ше боппаагу уаїше ргобіета 501уеа бу а (е8і їп м/рісП 18 ууає 
ипамаге ої Ше ехасі 50/ийоп соп5ігисіед и5іпє Ше (есппідцие 4езсгібед іп (|1|, а5 
арргоргіа/е, (гап5їегтед іо Ше са5е ої (гап5іепі Пеає сопдцспоп ргобіет. 

ІП Фе агіїсіе |2|, 18 15 соп5ідегед іп плоге декаї! Ше 5(гисіиге ої ре арргохітаге 
5оийоп, уПпеп Ше агеа 15 а гескапеіе. Пі Фі8 рарег, Ше сарабійцез Їог Ше са5е ОЇ гебіоп5 
соп5і5ШпЯ ОЇ гесіапеЇез аге іпуе5іеатей. Пі рагіїсиіаг, Ше ідеа ої їБі5 плеод 15 иц5ей іо 
соп5ігисі ап ехасі 50/ийоп ої ип5ісаду Преаї сопдисйоп ргобіепт, упіср Беіопе5 іо Ше 


сіавє СР" (р х(0, -юоо)) 
Тре ргобіет ої ип5ісаду реа сопдисййоп 15 соп5ідегед їп ап агеа Ша! ра5 Ше їогт 


ої сраппе! Баг58. 
Ап о оц зоГацоп 15 ве їп Фе Їогт ої: 


т-)| п-і 


(3 у,г) -УУс.9 нУУвид ук СА )Н, (5), 


ка (еп, кеть і (епі 
ууУреге (а у) є р, 120. 
Тре Типсйоп5 С,, (1) аге ТГоцпа Бу тіпіптіліпє Ше ехрге8зіоп / (п) изіпеє Фе 


тешфоад ої Їоса! роїепнаї: 


5 ЕВ 2) (2-х у, зн 


Трі8 Ісадз о а 5у5іет ої ппеаг ФНегеппа! едиайопя їп Ше уагіабіез С,, (1) 

По соп5ігисійпе Ше ехасі 50Їипоп Тог а (е5і са5е, Шеге пли5і Бе соп5ідегед Шаг Ше 
ехасі 50Гайоп Беопе58 іо Ше сраза ої Типсйоп5 Шаг Баме а сопійпиоц5 дегімайуєе їп 8 апа 
сопійпиоця 5есопа дегіуайуєзя и" и), . То до із, уе Пп Фе ехасі 50Їийоп ОЇ (е5і са5е 


іп Фе Рог ої 5 фНегепі Гогттиіаз їп 5 діїНегепі 5ибагеаз. 
Тре тешоа ргороведй їп Ше у/огК Тог соп5ігисіте, ехасі 50/ийопя ої (гап5іепі Пеаї 
сопдисйоп ргобіет5 їп агеа8, урісп сотрозеай ої гесіапеіез, аПоуму8 о ц5е Шезе 50ЇиОп5 


Бог сазе8 мбеге и(х, у, г) є С""! (рх (0, зюо)). 


Стаття надійшла до редакції 02.12.2011. 
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